Chromo-Electric flux tubes 
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The profiles of tfie cliromo-electric field generated by static quark-antiquark, QQ and three-quark, 
QQQ sources are calculated in Coulomb gauge. Using a variational ansatz for the ground state, we 
show that a flux tube- like structure emerges and combines to the "Y" -shape field profile for three 
static quarks. The properties of the chromo-electric field are, however, not expected to be the same 
as those of the full action density or the Wilson line and the differences are discussed. 
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I. INTRODUCTION 

An intuitive picture of quark-gluon dynamics emerges 
in the Coulomb gauge, V ■ A° = [1-3]. In this case 
QCD is represented as a many-body system of strongly 
interacting physical quarks, antiquarks and giuons. In 
particular the gluon degrees of freedom have only the two 
transverse polarizations and in the non-interacting limit 
reduce to the physical massless plane wave states. In 
the interacting theory gluonic states, just like any other 
colored objects, are expected to be non- propagating, i.e. 
confined on the hadronic scale. The non-propagating na- 
ture of colored states follows from the infrared enhanced 
dispersion relations which can be set up in the Coulomb 
gauge [3-6]. 

In the Coulomb gauge the A'^ component of the 4- 
vector potential results in an instantaneous interaction 
(potential) between color charges. Unlike QED, where 
the corresponding potential is a function only of the rel- 
ative distance between the electric charges, in QCD it is 
a functional of the transverse gluon components, A [1]. 
Thus the numerical value of the potential cannot be ob- 
tained without knowing the correct wave functional of the 
state and its dependence on the gluon coordinates. So 
in QCD the chromo-electric field is expected to be non- 
local and to depend on the global distribution of charges, 
which set up the gluon wave functional. 

Even though the exact solution to the general many- 
body problem is unavailable it is often possible to obtain 
good approximations if the dominant correlations can be 
identified. In Coulomb gauge QCD (in the Schrodinger 
field representation) the domain of the transverse gluon 
field, A is bounded and non- flat, and is referred to as 
the Gribov region. It is expected that the strong inter- 
action between static charges originates from the long- 
range modes near the boundary of the Gribov region, 
the so called Gribov horizon. For example it has been 
recently shown that center vortices, when transformed 
to the Coulomb gauge, indeed reside on the Gribov hori- 
zon [7]. 

The curvature of the Gribov region contributes to ma- 
trix elements via the functional measure determined by 
the determinant of the Faddeev-Popov operator. This de- 



terminant prevents analytical calculations of functional 
integrals, however it has been shown that its effect can 
be approximated by imposing appropriate boundary con- 
ditions on the gluon wave functional [9, 10]. This wave 
functional is in turn constrained by minimizing the ex- 
pectation value of the energy density which leads to a set 
of coupled self-consistent Dyson equations [4, 11]. Once 
the wave functional is determined it is possible to cal- 
culate the distribution of the chromo-electric field in the 
system. This is the main subject of this paper. 

In the following we study the chromo-electric field in 
the presence of the static quark-antiquark and three- 
quark systems, prototypes for a meson and a baryon 
respectively. Recent lattice computations indicate that 
the gluonic field near the static Q — Q state forms flux 
tubes. There are also indications that for the QQQ state 
the fields arrange in the so called "Y" -shape [12, 13], al- 
though some work supports the "A" -shape [14]. String- 
like behavior has been observed in the chromo-electric 
field in Ref. [15] and the "Y" -shape interaction advocated 
in Ref. [16]. A recent reevaluation of the center- vortex 
model also supports the "Y"-shape [17]. 

In the following section we summarize the relevant el- 
ements of the Coulomb gauge formalism and discuss the 
approximations used. This is followed by numerical re- 
sults and outlook of future studies. There is a funda- 
mental difference between lattice gauge flux tubes corre- 
sponding to the distribution for the action density and 
the chromo-electric field profiles. In the context of the 
potential energy of the sources, this difference was em- 
phasized in Zwanziger, Greensite and Olejnik [7, 18]. We 
discuss those in Section. IV. 



II. CHROMO-ELECTRIC COULOMB FIELD IN 
THE PRESENCE OF STATIC CHARGES 

A. The Coulomb gauge Hamiltonian 

The Yang-Mills Coulomb gauge Hamiltonian in the 
Schrodinger representation, H = -ff (11, A) is given by 

/'dx[n"(x).n«(x)+B'^(x).B''(x)]+Vc. (1) 
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The gluon field satisfies the Coulomb gauge condition, 
V • A°(x) = 0, for all color components a = 1 • • • — 
1. The conjugate momenta, n°(x) = — z£?/9A"(x) obey 
the canonical commutation relation, [11*'" (x), ^^'^(y)] = 
-i<5a(,(5^(Vx)<5(x - y), with ,5^(V) = 5^ - V^Vi/V^. 
The canonical momenta also correspond to the negative 
of the transverse component of the chromo-electric field, 
n''(x) = -E^(x), V • = 0. The chromo-magnetic 
field, B contains linear and quadratic terms in A. It will 
also be convenient to transform to the momentum space 
components of the fields by 



A"(k) = J rfxA"(x)e 



(2) 



and similarly for n"'(k). The Coulomb potential Vc may 
be expressed in terms of the longitudinal component of 
the chromo-electric field, 



Vc = \l dxE''(x)E"(x), 



with 



E«(x) = j 



dydz 



X - y 

47r|x — y|3 



1 - A 



(3) 



p'iz). (4) 



Here (1 — A) is the Faddeev-Popov (FP) operator which 
in the configuration-color space is determined by. 



where 



Aa5(p, q) = igfacb 



A^p - q) • q 



(6) 



/ are the SU (Nc) structure constants, and g is the bare 
coupling. In Eq. (4), p is the color charge density given 

by 



with the two terms representing the quark and the gluon 
contribution, respectively; the former is replaced by a c- 
number for static quarks. Without light flavors there is 
no other dependence on the quark degrees of freedom. 
The energy of the static QQ or QQQ systems measured 
with respect to the state with no sources is thus given 
by the Coulomb term and is determined by the expecta- 
tion value of the longitudinal component of the chromo- 
electric field. 

It is the dependence of the chromo-electric field and 
the Coulomb interaction on the static vector potential 
(through A) that produces the differences between QCD 
and QED. In QED the kernel in the bracket in Eq. (4) re- 
duces to [• • ■ ] S{y — z) and the Abelian expression for 
the electric field emerges. In QCD the chromo-electric 
field and the Coulomb potential are enhanced due to 



long-wavelength transverse gluon modes on the Gribov 
horizon where the FP operator vanishes. The combina- 
tion of two effects on the Gribov horizon: enhancement 
of (1 — A)^^ in the longitudinal electric field and van- 
ishing of the functional norm, which is proportional to 
det(l — A), leads to finite, albeit large, expectation val- 
ues of the static interaction between color charges. In 
Eq. (1) we have omitted the FP measure since, as men- 
tioned earlier in Ref. [9], its effect can be approximately 
accounted for by imposing specific boundary conditions 
on the ground state wave functional. 

Since the chromo-electric field depends on the distri- 
bution of the transverse vector potential it is necessary to 
know the wave functional of the system. A self-consistent 
variational ansatz can be chosen in a Gaussian form. 



"^[A] = exp 



dp 

(2^ 



■a;(p)A«(p) • A«(-p) . (8) 



The parameter uj{p) {p = |p|) is determined by minimiz- 
ing the expectation value of the energy density of the 
vacuum {i.e. without sources). The boundary condition 
referred to above corresponds to setting lj{0) = /i to be 
finite, which plays the role of Aqcd , *• e it controls the po- 
sition of the Landau pole. Minimizing the energy density 
of the vacuum leads to a set of coupled self- consistent in- 
tegral equations: one for lu, one for the expectation value 
of the inverse of the FP operator, d{p), 



{2wfS{p - cOSabdip) = 



I 



1- A 



ab 



|*)/(*|*), (9) 



and one for the expectation value of the square of the 
inverse of the FP operator, which appears in the matrix 
elements of Vc, 



{2nr6{k-q)6abf(p)d'{p) = 



p»(x) = V^(x)TV(x) + /„6cA^(x) • n^(x), (7) I 



1-A 



ab 



m/m. (10) 



J x,y 



The approximation / = 1 ignores the dispersion in the 
expectation value of the inverse of the FP operator, 



1-A 
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1-A 



(11) 



This approximation has been extensively used, e.g. in 
Refs. [2, 19]. The three Dyson equations were ana- 
lyzed in Ref. [4] where it was found that the solution 
of Lu can be well approximated by the simple function 
uj{p) = 6{fi—p)fj.+6{p—iJ,)p. The renormalization scale /i, 
being the only parameter in the theory, can constrained 
by the long range part of the Coulomb kernel (Vc) oc fd^. 
We will discuss this more in the subsection below. The 
low momentum, p < ji dependence of d{p) and of the 
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Coulomb potential Vc{p) = f{p)d{pf' is well approxi- 
mated by a power-law, 

dip) = din) [^y Jip) = /(m) {^^y (12) 

with a ^ 0.5 and {3 ^ 1. The exponents are boimded 
by 2a + /3 < 2 and the upper limit corresponds to the 
linearly rising confining potential. At large momentum, 
p » /i, as expected from asymptotic freedom, both d 
and / are proportional to l/log"'(p), with 7 = 0(1). 
Adding static sources does not modify the parameters 
of the vacuum gluon distribution, e.g. u){p). This is 
because the vacuum energy is an extensive quantity while 
sources contribute a finite amount to the total energy. 
Thus we can use the three functions u), d and / calculated 
in the absence of the sources to compute the expectation 
value of the chromo-electric field in the presence of static 
sources. The ansantz state obtained by applying quark 
sources to the variational vacuum of Eq. (8) does not, 
however optimize the state with sources. 



B. The field lines in the QQ and QQQ systems 

For a quark and an antiquark at positions = R/2 = 
Rz/2 and Xg = — R/2 = —Rz/2, respectively, and the 
gluon field distributed according to ^'[A], the expecta- 
tion value of the square of the magnitude of the chromo- 
electric field measured at position x is given by 



^^'('^'^)> = 7S^ ^ ^ ±fdy,dy, 

^ ' zi=±R/2 Z2 = ±R/2 

(x-yi) ■ (x-y2) 



|x - yi|^|x - y2|^ 



-E(zi,yi;z2,y2), (13) 



where the +(— ) sign is for the zi = (7^)z2 contributions, 
and 



-E'(zi,yi;z2,y2) = 
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g 




1-A 




1-A 










y2,Z2 



I*) 



(14) 



The color factors leading to Cp can be extracted from 
the expectation value in Eq. (14) (the ground state ex- 
pectation value of the inverse of two FP operators is an 
identity in the adjoint representation). In the Abelian 
limit, Ei^i ■ • • y2) 5{yi - zi)(5(y2 - Z2) and Eq. (13) 
gives the dipole field distribution, {E?)qed- One should 
note that Eq. (13) contains the two self energies. These 
self energies are necessary to produce the correct asymp- 
totic behavior at x >> i? for charge- neutral systems, (in 
QED and QCD) i.e has to fall-off at least as 1/x^ at 
large distances from the sources. 

The infrared, |x| ~ |R| >> 1//^ enhancement in QCD 
arises from the expectation value of the inverse of the FP 
operator. If (E^(x, R)) is integrated over x one obtains 



the expectation value of the Coulomb energy of the QQ 
source. The mutual interaction energy is given by. 



Vb(R) = ^ j rfx(E2(x,R)) 



9 



-Cf 



1 - A V V2 y 1 - A 

dp dHp)f{p) 



ip R 



(27r)3 p2 
and the net self-energy contribution is. 



1 

V2 



1-A 
dp d\p)f{p) 



9 



1-A 



(27r)3 



p- 



(15) 
(16) 



In lattice simulations it has been shown [20] that the 
Coulomb energy and the phenomenological static QQ 
potential obtained from the Wilson loop are different. 
In particular it was found that the Coulomb poten- 
tial string tension is about three times larger than the 
phenomenological string tension. This is in agreement 
with the "no confinement without Coulomb Confine- 
ment" scenario discussed by Zwanziger [18]. It is sim- 
ple to understand the origin of the difference. Even if 
|\1/[A]) were the true vacuum state (without sources) of 
the Coulomb gauge QCD Hamiltonian (here we approx- 
imate it by a variational ansatz) the state \QQ,R) = 
(5(R/2)(5(-R/2)|*[A]) is no longer an eigenstate. For 
example Vc acting on \QQ,R) excites any number of 
gluons and couples them to the quark sources. The 
Coulomb energy was defined as the expectation value, 
Vc in \QQ, R) minus the vacuum energy and it is there- 
fore different from the phenomenological static potential 
energy which corresponds to the total energy (measured 
with respect to the vacuum) of the true eigenstate of the 
Hamiltonian with a QQ pair. If one defines [7] 

G{R,T) = {QQ,R\e-("-^»^^\QQ,R) 

= I {QQ, R, n\QQ, i?) |2e-(^"-^°)^, (17) 

n 

then the Coulomb potential on the lattice can be calcu- 
lated from 



Vc{R) = hm -— log(G(i?,T)), (18) 
and the phenomenological potential from 



y(i?)= lim - log(G(i?,T)). 

T=oo dl 



(19) 



Thus one should be comparing Vc{R) in Eq. (15) to the 
lattice Coulomb potential and not to the phenomenolog- 
ical potential obtained from the Wilson loop. Finally, 
one could try to optimize the state with sources, e.g. 
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by adding gluonic components. In this case terms in the 
Hamiltonian beyond the Coulomb term would contribute 
to the energy of the system and one could compare with 
the true (Wilson loop) static energy. In our previous 
studies, where we extracted numerical values for /x and 
the critical exponents a,/?, 7 (c/. Eq. (12)) we have in- 
stead compared Vc to the phcnomcnological, Wilson po- 
tential [4] . In what follows we will use the larger value of 
the string tension, to be in agreement with Rcf. [7] . 

If the two exponents a and /3, which determine the 
infrared behavior of d{p) and f{p) respectively, satisfy 
2a + P > 2, then the self energy in Eq. (16) is divergent 
and so is the rhs of Eq. (15). This reflects the long-range 
behavior of the effective confining potential generated by 
self-interactions between the gluons that make up the 
Coulomb operator. For the colorless QQ system the total 
energy which is the sum of Vc and S, is finite as it should 
be. For a colored system, e.g. a quark-quark source, the 
sign of Vc changes, there is no cancellation between the 
infrared singularities, and in the confined phase the sys- 
tem would be un-physical with infinite energy. The inte- 
gral determining the self energy also becomes divergent 



in the UV, since for p ^ 00 the product (P{p)f{p) only 
falls-off logarithmically. Modulo these logarithmic cor- 
rections this UV divergence is the same as in the Abelian 
theory and can be removed by renormalizing the quark 

charge. 

It follows from translational invariance of the matrix 
element in Eq. (14), that E depends only on the relative 
coordinates, zi — yi and Z2 — y2. We therefore introduce 
the momentum space representation, 



-E(zi,yi;z2,y2) 



!^_^^gip-(zi-yi)-iq-(z2-y2) 



(27r)3 (27r)3 
xd{p)d{q)E{p;q) 



(20) 



and define Fl(1) = E{1 + L/2;l - L/2) with 1 = (p 

q)/2 and L = p — q. The Dyson equation for F can be 
derived in the rainbow-ladder approximation which, as 
shown in Ref. [4, 5], sums up the dominant infrared and 
ultra-violet contributions to the expectation value of the 
inverse of two FP operators, 



J 



FL(l) = l + iVe 



dk [(k - L/2)5T(k + I)(k + L/2)] d{k - L/2) d{k + L/2) 



(27r)3 



2a;(k -|- 1) 



(k-L/2)2 (k + L/2)^ 



^L(k), 



(21) 



r 



It follows from Eq. (20) that L and 1 are conjugate to the 
center of mass, R = [(zi — yi) -|- (z2 — y2)]/2 and the 
relative, r = [(zi — yi) — (z2 — y2)] coordinate respec- 
tively. The Dyson equation for Fj^ is UV divergent if for 
p/fj, » 1, and d{p) > log^/^(p^). This divergence can be 
removed by the Coulomb operator renormalization con- 
stant. The renormalized equation is obtained from the 
once-subtracted equation Fl(1) — Flo(Io). For example, 
if the subtraction is chosen at |lo| = /i and Lq = 0, 
the renormalized coupling -fo(M) can be fixed from the 
Coulomb potential. After integrating Eq. (13) (over x) 



one obtains (5(yi — y2) multiplying E{zi, ■ ■ ■ , y2). There- 
fore, it follows from Eq. (20) that Vc(R) is determined by 
Fo(l) and Fo{l) = Fo{l) = f(l) with / defined in Eq. (10). 

In Eq. (21) L is a parameter, i.e. the Dyson equa- 
tion does not involve self-consistency in L. We have just 
shown that as L ^ 0, ^l(I) has a finite limit: it is given 
by /. For large L — |L| {L/iJ » 1), due to asymp- 
totic freedom, Fl is expected to vanish logarithmically, 
Fl d'^{L) oc l/log(L^). We do not attempt here to 
solve Eq. (21), instead we use a simple interpolation for- 
mula between the L = and L ^ 00 limits, 



J 



FL(l) = /(l)e(M- 



p + q 



9{n - p)e{fi - 9) + [1 - e{p)e{q)] , (22) 



r 



i.e. in the term in the bracket we ignore the short dis- 
tance logarithmic corrections. It is easy to show that if 
logarithmic corrections are ignored then the short-range, 
p,q > fj, contribution to the energy density is the same 
as in the Abelian case. Since we are mainly interested 



in the long range behavior of the chromo-electric field, 
in the following we shall ignore contributions from the 
region p,q > ^ all together. In the long-range approxi- 
mation, x,R» l//i the expectation value of is then 
given by. 
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(E2(x,R)) = 



(4 



^ V p'^^ f 

:7r)2 4- J 



drhir)-. 



-x-r/2| |zj-x + r/2| 



+ ""/^ 4(zi - X - r/2)di(z,- - X + r/2). (23) 



r 



eg-^ = 1 for i = j and -1 for i ^ j, Zi,(2) = (-)R/2 
2 /"^ 



pdpji{rp)d{p), 



is the derivative of w.r.t. r, 

1 /•A' 



^ ^ dpp^ f{p)jo{pr), 



(24) 



(25) 



and joiii ai"© Bessel's functions. We note that the ex- 



pression in Eq. (23) is not necessarily positive. In the 
limit f{p) = 1, the matrix element of the square of the 
inverse of the FP operator is approximated by the square 
of matrix elements (cf. Eq. (11)) and (E^) becomes pos- 
itive. 

The expression for (E^) for the three quark sys- 
tem is derived by taking the expectation value of 
the Coulomb operator, Vc in a color-singlet state 
ey7s(3i(zi)(5j(z2)(3fc(z3)|*[A]), which gives 



J 



^ ' ij=l 



X — r/2 Zj — x -|- r/2 , 



x-r/2| |zj-x-|-r/2| 



d'i(zi - X - r/2)d'i(zj - X + r/2) (26) 



where ^^'^'^ = i if ^ = j and ^fj'^'^ = -l/2iii^j. We 
note that the energy density for the QQQ system comes 
from two-body correlations between the QQ pairs. 



cQQQ 



III. NUMERICAL RESULTS 

We first consider the simple approximation to the ex- 
pectation value of the Coulomb kernel of Eq. (11) in 
which f{p) = l. If one wishes to have the confining 
potential grow linearly at large distances then it is nec- 
essary to set a = 1, i.e. d{p) oc ^/p for p/fj, < 1. In 
this case, assuming that the long-range behavior of the 
potential is of the form Vc{r) = bcr, we obtain from 
Eq. (15), 



bc = CFd^{^i)^iV{8^^). 



(27) 



Wo use the Coulomb string tension be = 0.6 GeV^. For 
the QQ system the long-range contribution to the electric 
fields is then given by. 



(E2(x,R)) = 



2bc 



(R/2-x) 



|R/2-x|2 



(l-jo(/x|R/2-x|)) 



+ (x ^ -x) 



(28) 



In Fig. 1 we show the Coulomb energy density as a 
function of position on the QQ axis, x = a;R, for R = 
|R| = 10 fm. The small oscillations come from the sharp- 
cutoff introduced by the ^-functions in Eq. (22) which 




FIG. 1: R^(E^{x)} in units of 2bc/n^{hcf as a function of 
the distance x along the QQ axis. We employ the f{p) = 1 
approximation. The quark and the antiquark are located at 
R/2 = 5 fm and —R/2 = —5 fm respectively. The renormal- 
ization scale /x = 1.1 GeV is calculated from Eq. (27) using 
d(At) = 3.5 from Ref. [4]. The dashed line is the contribution 
from the two self energies, the dash-dotted line represents 
mutual interactions and the solid line is the total. 



produces the Bessel's functions in Eq. (28). For a smooth 
cutoff, e.g. with ^(/u — p) — > exp{—p/ijb) in Eq. (28) one 
should replace 1 — jo(/i|R/2 — x|) by 1 — arctan(/i|R/2 — 
x|)//z|R/2 — x|. The cut-off is also responsible for the 
rapid variations near the quark positions, x = ±i?/2. 

We note that for large separations between the quarks, 
R » 1//X and x « R, the Coulomb energy density 
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FIG. 2: i?^(E^(a:)> in units of 2bc/Tt^{hc)^ as a function of 
the distance x transverse to the QQ axis. The units and the 
setting are as in Fig. 1. 



behaves as expected from dimensional analysis, 

326c 



(E^(x,i?/i oo)) 



(29) 



which is consistent with linear confinement, i.e. if 
(E2(x. oo)) is integrated over x in the region 

|x| < i? on obtains Vc{R) oc R. 
At large distances a;>>ii>>l/;uwe obtain 



(E^dxl/i? ^oo,Rfx^ oo)) 



2bcR' 



(30) 



If there were a finite correlation length one would expect 
(E^dxj/i? — »• 00, R^ — > oo)) to fall-off exponentially with 
|x| [8] and not as a power-law. The power-law behavior 
obtained in Eq. (30) is again related to the difference 
between the \QQ, R) state used here, which is built by 
adding quark sources to the vacuum and the true ground 
state of the QQ system as discussed in Sec. IIA. In other 
words the profile of the chromo-clcctric field distribution 
for such a state is not expected to agree with the pro- 
file of the fiux-tube or action density. To illustrate this 
difference, in Fig. 2 we plot the energy density as a func- 
tion of the magnitude of the distance transverse to the 
QQ axis, x± = |x_l|, R • x = R • xx =0. 

Finally, in Fig. 3, we show the contour plot of the en- 
ergy density as a function of the position in the xz plane 
with quark and antiquark on the z axis at R/2 and —R/2 
respectively. 

It is clear from Figs. 2 and 3 that a fiux tube like struc- 
ture emerges and from Eq. (29) that it has the correct 
scaling as a function of the QQ separation but, as dis- 
cussed above it does not have a finite correlation length 
(large x behavior). 

The field distribution for the QQQ system in the 
/l(p) = 1 approximation is equal to the sum of three 
terms each representing a contribution from a QQ pair. 
We place each of the three quarks in a corner of an equi- 
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FIG. 3: R'^{E'^{x)) as a function of position in the xz plane. 
The units and the same setting as in Fig. 1. 



lateral triangle, Zj 



Cf 



(E2(x,R,)) = ^[(Dr-D2) 



+ (Di-D3)2 + (D2-D3)2], (31) 



where 



D, 



■r/2 



|zj - x + r/2| 



r/2). (32) 



The contour plot of of energy density in this case is 
shown in Fig. 4. Even though the field originates from the 
two-particle correlations the net field seems to form into 
a "Y" -shape structure. This structure has also recently 
been seen to emerge in Euclidean lattice simulations. 

Finally, to study the effects of /l(p), in Fig. 5 we show 
the predictions for the QQ field distribution given by 
Eq. (15) where we use d{p) and f{p) in the form given 
by Eq. (12) with a = 1/2 and (3—1 and normalized 
such that V{R) — * bR at large distances. Furthermore, 
to remove the oscillations introduced by the momentum 
space cutoff, wc now cut the small x region in coordinate 
space, by i) extending the upper limits of integration in 
Eqs. (24) and (25) to infinity and ii) cutting off the po- 
sition space functions at short distances. 



d'^{r) = — sin(^a/2)r(2 - a)0{rn - 1)-^ 

TTQ (a**") 



Mr) = -^sin(7r/3/2)r(2 - /3)0(r/z - 1) 



27r' 



(/xr) 



3-/3" 



(33) 



(34) 



Comparing Fig. 3 and Fig. 5 we observe a narrowing 
of the flux tube. This is to be expected as the action 
of f{p) is to introduce additional gluonic correlations. 
That said, there is no major qualitative change in the 
field distribution. 
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FIG. 4: R'^{E'^{x)) as a function of position in the xz plane. 
The units and the same setting as in Fig. 1. The upper panel 
show the total field distribution and the lower the distribution 
from mutual interaction (no self-energies) only. 
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FIG. 5: R^{E'^{x)) for QQ from Eq. (23) with a = 1/2 and 
(3=1. The units and the same setting as in Fig. 1, except 
that the contribution from f(p) has been included. 



IV. SUMMARY 

We have calculated the distribution of the longitudi- 
nal chromo-electric field in the presence of static QQ and 



QQQ sources using a variational model for the ground 
state wave functional. Despite this wave functional hav- 
ing no string-like correlations a flux tube like picture does 
emerge. In particular the on-axis energy density of the 
QQ system behaves as bc/B? for large inter-quark sep- 
aration, R and the field falls off like bcR^/x'^ at large 
distances from the center of mass of the QQ system, x. 
This is weaker than in the Abelian case (~ /x^) and 
implies that moments of the average transverse spread of 
the tube, defined as proportional to {\x±_\^'E?'{z,x±_)), are 
finite for n < 2 only. Thus there is no finite correlation 
length for the longitudinal componc;nt of the chromo- 
electric field, as expected for the state which does not 
take into account screening of the Coulomb line by the 
transverse gluons (flux tube) . This also leads to large Van 
der Waals forces, which is bothersome, but it is consistent 
with the scenario of "no confinement without Coulomb 
confinement" of Zwanzigcr. The Coulomb potential leads 
to a variational (stronger) upper bound to the true con- 
fining interaction. 

Similar behavior at large distances is also true for the 
three quark sources, except that here we find the emer- 
gence of the "Y"-shape junction. This is consistent with 
lattice simulations, but is remarkable in our case as it 
arises from two-body forces. It will be interesting to ex- 
amine field distributions which include transverse field 
excitations. In that case the only lattice results avail- 
able are for the potential, not for the field distributions. 
Finally we note that, since the mean field calculation 
provides a variational upper bound, the long range be- 
havior of the field distribution falls-off more slowly than 
expected for the Van der Waals force. Certainly as the 
complete string develops this is expected to disappear 
and it would be interesting to build a string-like model 
for the ansatz ground state to verify this assertion. 
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